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■ In this paper we report exact results on thermostatic properties of the classi- 

Q\ . cal n-vector model on the harmonic chain. This system is characterized by the 

Hamiltonian 

H:=#o(M,{?}) {£}) (1) 
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Here Hq is the Hamiltonian for the nearest-neighbor coupled harmonic chain which 
consists of a one- dimensional lattice £Z with lattice constant £ > and a set of iV + 1 
point particles of mass m > distributed along the Euclidean line IR at positions 
j£ + qj, j = 1,2, N + 1. The momentum of the j-th particle is denoted by pj 
and the spring constant of this chain is itluJq. We assume now that each particle 
carries a set of internal rotational degrees of freedom which we collectively represent 
by a classical spin, that is, by an n-component Euclidean unit vector Sj G R n . The 
Hamiltonian Hi then models the simplest rotational invariant interaction between 
the spins of two nearest neighboured particles. The interaction strength between 
two spins is described by the real- valued even function W : x i— > W(x) and depends 
on the actual interparticle distance as indicated in (3). In accordance with the 
harmonic approximation it is sufficient to consider only the first two terms in a 
Taylor expansion of W, 

W(£ + q j+1 - qj ) w J + (q j+1 - qj )r], (4) 

where J := W(£) and rj := W'{£). 

Within this approximation it is possible to decouple the vibrational and rota- 
tional degrees of freedom by introducing shifted particle positions (cf. |], 131 f° r 

71=1) 

i-i 

x x \=qx, x 3 :=q 3 - j S r • S r+ i , j = 2, 3, . . . , N + 1. (5) 

muJ r= i 
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The result can be cast into the form 

H « H ({p},{x}) + # spin ({£}) (6) 
where we have introduced the pure spin-chain Hamiltonian 

JV 

H spm ({S}) := - £ (j$ ■ + ■ 5 J+1 ) 2 ) (7) 
i=i 

with X := rf/(2mujQ). Since the thermal properties of the harmonic chain are well 
known (see, for example, 0), we will consider only those of H spin . 

Special cases of the Hamiltonian (7) have already been discussed in the litera- 
ture. For K — it corresponds to Stanley's n- vector model in one dimension ||, [5]. 
For J = and n G {2,3} a discussion is due to Vuillermot and Romerio ||. As 
for the Ising case (n = 1), we remark that the biquadratic term in (7) lowers the 
specific free energy of the Ising chain simply by the constant K. However, even for 
n — 1 this term is responsible for magnetostrictive effects of the full system (6), as 
discussed by Mattis and Schultz ■ As an aside we mention that exactly known 

|7| ground-state properties of the quantum version of (7) for n = 3, Sj = 2h 2 l, 
J < and K/\J\ = — 1/3 are discussed in relation with Haldane's conjecture. 



The basic thermostatic properties of the classical spin-chain Hamiltonian (7) can 
be obtained from the free energy per spin in the macroscopic limit iV — ► oo: 

where the canonical partition function at temperature l/(3k (k: Boltzmann's con- 
stant) for the finite chain and n > 2 may be defined by the (N+l)(n—l) -dimensional 
integral 

Z(f3) := JdSi---J dSW+i exp {-/3# spin ({£})} . (9) 

For convenience we are using open boundary conditions. Furthermore, each of the 
above dS* stands for the usual surface measure on the (n— l)-dimensional unit sphere 
in R n . We assume this measure to be normalized in the sense that / dS = 1 and 
recall its invariance under rotations 

J dSf(S) = J dSf(gS). (10) 

This relation is valid for any integrable complex- valued function / and all orthogonal 
n x n matrices g G SO(n). Of course, for n = 1 the integration / dS stands for the 
summation ^J2s=±i- 

For the evaluation of the partition function (9) we note that the Hamiltonian (7) 
can be rewritten as H spm (^{S}^j = J2jLi V0j, Sj+i) where we have introduced the 
spin-pair interaction energy 

V{S,S') := -JS-S' -K(S-S') 2 , (11) 
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which is 50(7?,) -invariant and exchange- invariant: 

V(gS, gS') = V{S, S') = V(S', S), for all g G SO{n). (12) 



These properties can be used to rewrite the Hamiltonian (7) as H spin y{S}^ 
J2jLi V(So, gjSj+i) where So is an arbitrary but fixed unit vector and the n x n 
matrices gj are defined by gjSj := Sq. With the rotational invariance (10) the 
partition function (9) can be reduced to a single dS-integration according to 



Z(p) = J dS, j dS 2 e- f3V ^ s ^ ■ ■■JdS N+1 e~l 3v ^ N +^ = \ N (f3) (13) 

where 

\({3) := J dS ex V {-(3V(S , S)}. (14) 

Hence, the specific free energy (8) is given by 

F((3) = -(l/(3)\n\((3). (15) 

What remains to be done is the integration (14). Choosing as the fixed vector 
So the unit vector pointing towards the northpole, So = (0, . . . , 0, 1), the function 
V(Sq,S), and therefore also exp{— PV(So, S)}, depends only on the polar angle 
9 if S is parameterized in the usual (hyper-) spherical polar coordinates, because 
then So ■ S = cos 9. Hence, the expression (14) can immediately be reduced to the 
following one-dimensional integral (t := cos 9): 

A03) = -P^tt / dte^+- 2 )(l - t 2 )^. (16) 
For K = this integral can be expressed in terms of modified Bessel functions 



n — 2 
2 



A(/3) = T(n/2) (Jj j l^J), (17) 

the well-known result for Stanley's n- vector chain ||. For K ^ an expansion 
in powers of (3 J allows for an integration in terms of the confluent hypergeo metric 
function ||: 

This series can be summed M in terms of a generalized hypergeometric function of 



two variables 10, 11 



A(/3) = exp {-#} * 2 (i; % ±; (3K, . (19) 

Unfortunately, not much is known about this generalized hypergeometric function. 
However, for the cases n — 1, 2 and 3, we can express A(/3) somewhat more explicitly 
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in terms of hyperbolic, modified Bessel and confluent hypergeometric functions, 
respectively: 



A(/3) 
A(/3) 
HP) 



e pK cosh(/3 J) 



b : " J ° 2/4K [(i + *) ^ 



for 77, = 1, 



for n 



for n 



(20) 



For n = lwe have used formula 7.2.4.91 of Ref. |L0[ leading to the expected result. 
For n = 2, in essence, we have expanded the integral (16) in powers of (3K. For 
n = 3 the integral (16) obviously is reducible to the sum of two error functions 
with complex argument, which in turn are expressable in terms of confluent hyper- 
geometric functions. Eqs. (19) and (20) in combination with (15) summarize the 
main results we wish to report here. Appropriate derivatives of A(/3) with respect 
to (3 lead to the basic thermostatic quantities. For example, the specific heat c(/3) 
is given as 



rj2 



d(3 



"A"(/3) 






{ A(/8) J _ 



(21) 



Figure 1 displays this function for n = 3 as a surface over the l//3\ J\-K/\ J|-plane. 
We note that its zero-temperature value c(oo) = k is in agreement with the classical 
equipartition theorem. For low temperatures and K > the specific heat c(f3) 
increases linearly: 

c((3) , K 



1 + 



+ 



k " ' (K+|J|/2)2/3 
For high temperatures it vanishes as the inverse square of the temperature: 



k 



J 2 + — K 2 
15 



P 2 



+ 



(22) 



(23) 



For < K/\J\ ^ 14 the specific heat attains a maximum value in the temperature 
range < 1/(3\J\ ^ 0.9. This maximum splits for K/\J\ 15 into two maxima. 
The global maximum of c(/3) remains near 1//3\J\ ^ 1 for large K values. 



Finally, we remark that from the free energy (15) in combination with (16) one 
can obtain |E| further interesting properties of the system (6). Examples are the 
two-spin correlation function 



(24) 
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here (•) denotes the canonical equilibrium expectation value with respect to Hamil- 
tonian (6), the zero- field susceptibility 

xM ■= P (i + 2 E(4 • 4*)) = P \ + £ F F {py ( 25 ) 

and the mean lattice constant 

aifi) := t + (q j+1 - qj) = £ - -\ |- F{0). (26) 

A more detailed discussion (including applications and further thermostatic prop- 
erties) of the system characterized by the Hamiltonian (6) will be given elsewhere 

El- 
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Figure Caption 



Figure 1: The specific heat (21) for n = 3 as a function of the dimensionless tem- 
perature kT/\J\ := 1/P\J\ and the parameter K/\J\. 
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